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In a series of papers Boyanovsky et al. have studied the evolution of an 



' inflaton with a negative mass squared and a quartic self coupling using the 

Closed Time Path (CTP) formalism relevant for out-of-equilibrium dynamics. 
' In this paper we comment on various aspects of these works. We first compare 

' their approach to alternate approaches to study inflaton dynamics and point 

Q_i' out that the use of the CTP formalism gives the same results as standard field 

theory in the Hartree and leading order large N approximations. We then 



rederive using the WKB approximation the large momentum mode functions 



^ ' of the inflaton needed for renormalisation and point out some differences with 



the previously obtained results. We also argue that the WKB approximation 
is valid only for large k/a and not for large k as apparently assumed in the 
above mentioned works. We comment on the renormalisation prescription 
adopted in these works and finally discuss how it differs from another more 
commonly used prescription. 
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I. INTRODUCTION 



In the new inflationary scenario the expansion of the universe is driven by the nearly 
constant potential energy of a slowly rolling scalar field. In Refs. [1,2], Boyanovsky et al. 
have discussed in great detail the evolution of an inflaton with a negative mass squared and 
a quartic self coupling using the Schwinger-Keldysh Closed Time Path (CTP) formalism 
[3,4] which is suited for non-equilibrium or dynamical situations. The advantage of the 
CTP formalism is that it allows one to obtain (in\in) matrix elements of field operators (as 
opposed to (out\in) matrix elements of standard field theory) which are relevant in non- 
equilibrium or dynamical situations [5]. 1 In an earlier work [7] it had been shown (in 
the context of Minkowski spacetime) that if unstable modes grow large then a perturbative 
treatment within the CTP formalism will be insufficient and one needs to invoke a non- 
perturbative treatment. Therefore in Ref. [1,2] both the Hartree and the large N non- 
perturbative approximations are used to obtain the equation of motion for the inflaton. 

Renormalisation was achieved in Refs. [1,2] by subtracting off the divergent ultraviolet 
behaviour based on the renormalisation scheme employed in Ref. [6]. The ultraviolet be- 
haviour was obtained by analytically solving the equations of motion for the large momentum 
fluctuations of the inflaton field after applying the Hartree or large iV approximations. These 
approximations render the equations of motion linear and thus amenable to techniques such 
as the WKB method. 

In this work we first restate the equations of motion and the expression for the inflaton 
fluctuations obtained in Refs. [1,2]. We then compare this with earlier approaches to studying 
inflaton dynamics. We next rederive the expressions for the large k mode functions of the 



1 Alternatively, in Ref. [6] the authors use a Gaussian ansatz for the density matrix, which is 
evolved in time with a time dependent Hamiltonian using the von Neumann equation, and this is 
used instead of the CTP formalism to obtain the necessary expectation values for the dynamical 
system. 
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inflaton field needed for renormalisation, after introducing more generalised initial conditions 
with finite temperature corrections to the initial mass of the inflaton. Since some of our 
results differ from those in Refs. [1,2] and related works, and since the calculations are 
involved, we have provided all the details and pointed out where we differ from previous 
works. We have also commented upon certain aspects of the renormalisation procedure. We 
then compare the renormalisation scheme adopted in Refs. [1,2] with schemes adopted by 
other authors. Lastly we provide some numerical results. 

To study the evolution of the inflaton in its potential we choose initial conditions for 
the field modes consistent with those specified in Ref. [2]. We work in the context of the 

— * 

large N approximation and introduce an iV-plet $ to represent the inflaton. As in Ref. 
[1] we work in the context of a quenched approximation in which the initial state of the 
universe at the onset of inflation at t — corresponds to a state of positive mass $ particles 
in thermodynamic equilibrium at a temperature Tj > T c , the critical temperature, while 
the subsequent evolution is for a universe at T = in which the mass term for <3> becomes 
negative. We assume that the distribution of <£> particles freezes out soon after the onset of 
inflation at a temperature close to Tj. 

Our major conclusions are the following: 

• The equations of motion obtained in Refs. [1,2] afer invoking the CTP formalism are 
the same as those obtained using standard field theory. This is because of the Hartree 
and lowest order large N approximations adopted in Refs. [1,2]. 

• The WKB solutions obtained in Refs. [1,2] for large k (i.e. k 3> m R , where m R is the 
renormalised inflaton mass parameter) are actually valid for large k/a (i.e. kja~^> m R ) . 
Since the renormalisation scheme of Refs. [1,2] involves subtracting the contribution 
of modes with k > m R the counterterms used in Refs. [1,2] do not match the bare 
quantities for m R a > k > m R . 

• The dominant modes relevant for inflaton evolution in new inflation where the inflaton 
has a negative mass squared are the low momentum modes. The counterterms used 
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in Refs. [1,2] do not have any significant effect on the contribution of these modes and 
hence, despite our above criticism of the counterterms, one gets the same results as in 
the more standard prescription for removing ultraviolet divergences for the inflaton. 

• In the more standard prescription to deal with ultraviolet divergences for the inflaton 
one subtracts the contribution of modes with k/a > H. In Refs. [1,2] the intent of 
the renormalisation scheme appears to be to subtract the contribution of modes with 
k > m R . While for the case of an inflaton with a negative mass squared it turns out 
that both schemes give the same results, for an arbitrary field the renormalisation 
scheme of Refs. [1,2] should be used with care. We illustrate this with an example of 
a massless field in de Sitter space. 

II. EQUATIONS OF MOTION AND INITIAL CONDITIONS FOR INFLATON 
MODES, AND TOTAL FLUCTUATIONS 

The Lagrangian density that we consider in the context of a spatially flat Robertson- 
Walker universe is given by, 



£ = a 3 (t) 



(1) 



$ is an iV-plet and M 2 = — m 2 + £ 1Z. 1Z is the Ricci scalar and is given by 



\a(t) a*(t)) 

The scale factor can be obtained dynamically using the Einstein equation. (In the numerical 
solutions in Ref. [2] H was obtained dynamically while it was taken as a constant in the 
simulations of Ref. [1].) 

Below we sketch the large N approximation as discussed in greater detail in Refs. [8,1]. 
We then present the equations of motion and discuss the initial conditions. In this section 
we also discuss other approaches to study inflaton dynamics. 



Non-perturbative treatment of non-equilibrium dynamics can be studied by approxima- 
tions such as the Hartree approximation and the large N expansion [1]. These approxima- 
tions are similar at leading order in 1/N. Here we shall work with the large N approximation. 
In the large iV approximation the Lagrangian can be written in terms of X/N as above and 
in the spirit of Ref. [1] it is assumed below that A is small enough that the large N approx- 
imation is valid even when N = 0(1). (We shall choose A = 1CT 12 later.) Following Ref. 
[8,1], only leading order terms in the large N approximation are retained below. 

To facilitate the use of the large N approximation to study the inflationary phase transi- 
tion the inflaton is written as 

$(x,t) = (a(x,t),7f(x,t)), 
where n is an N — 1-plet, and we let 

a(x,t) = cr (£,*) +p(x,t) ; (cr(x,t)} = a (x,t) ; (p(x,t))=0. (4) 
To apply the large N approximation <r and tt are defined by 

N-l 

a (x, t) = y/N(/>(t) , 7?(£ , t) = ^(f , t) (1~V~~I) , (5) 

mimics the mean of the inflaton field while ip and p correspond to fluctuations. It is 
implicitly assumed here that the mean field is a function of time only. 
The mode functions £4(t) of the inflaton are defined as 

with [ak, ajj/] = (27r) 3 5(k — k'). The leading order in 1/N equations of motion for and 
Uk(t) then are 

0(t) + 3^(0 + M^(t) + ^(t) + ^{m 2 {t)) = 0, (7) 

U k (t) = 0, (8) 
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and the effective frequency u>k(t) is given by 



<4® 



+ M 2 (t) , 



(9) 



where 



M 2 (t) = M| + -0 2 (t) + -(^(t)), 



(10) 



and 



d 3 k \U k (t)\ 2 
(2tt) 3 2 



coth[W fc0 /(2T,)]. 



(11) 



By (ip 2 (t)} we mean {^ 2 (x,t))(t). By translational invariance (ijj 2 (x,t))(t) depends only on 
time. 

The coth function reflects the initial thermal particle distribution and it is assumed that 
ip goes out of equilibrium at a temperature close to Tj. We obtain Tj by specifying that 
inflation starts when the energy density in radiation falls to one-tenth of the constant term 
in the potential energy. Assuming the number of relativistic degrees of freedom to be 100, 
we get Ti ~ 200m R . Wko is related to the initial frequencies of the modes and is defined 
later. 

At leading order in 1/N, p does not enter in the equations of motion for <fi and Uk- 
Furthermore the contribution of p to the energy density is smaller by a power of iV compared 
to that of the other fields. Hence it is not included in the analysis. 

The above equations of motion are similar to those in earlier works that do not invoke 
the CTP formalism, such as Ref. [9]. Modifying Eqs. (2.13) and (2.14) of Ref. [9] for our 
potential (and our definition of the Fourier modes) gives 




(12) 



where 




(13) 
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The similarity in our equations and those of Ref. [9] 2 , though the latter are obtained without 
invoking the CTP formalism, is because we are working to lowest order in 1/N, and at this 
order both the (out\in) and the (in\in) formalisms give the same equations of motion. This 
may be verified by looking at the field equations obtained using standard field theory and 
the CTP formalism in Ref. [5]. Eqs. (4.27)-(4.34) of Ref. [5] give the equations of motion 
in the CTP formalism including terms of order 1/N for the mean field and the momentum 
modes and these reduce to equations obtained using standard field theory, namely, Eqs. 
(2.20) and (2.21) of Ref. [5] for the mean field and Eq. (12) above for the momentum 
modes, if we ignore terms of order 1/N. 3 Note that at lowest order in 1/N the equal time 
Wightman functions i(vp 2 (t)) obtained using the CTP formalism (also referred to as G > ,G < ) 
and the equal time Feynman propogator obtained in standard field theory that appear in 
the equations of motion are the same (see Sec. II of Ref. [10]). 

Another approach to studying inflaton dynamics is provided in Ref. [11] in which the 
mean field (<&) is taken to be and the dynamics is studied by evolving ($ 2 ) by defining 
an effective action which is a function of ($ 2 ) rather than of (<&), similar to the Cornwall- 
Jackiw-Tomboulis method [12] in which one studies the two particle irreducible effective 
action. In Ref. [11] the authors do not invoke the CTP formalism. In Ref. [13] the CJT 
method is combined with the CTP formalism to study the evolution of a A0 4 scalar field 
with a tachyonic mass which is suddenly brought into contact with a heat bath at zero 
temperature (but in Minkowski spacetime). The CJT method directly gives an equation of 
motion for the two-point function. In contrast, above, as in Refs. [1,2,5], we obtain the time- 

2 The factor of 3 with A above is because the cactus approximation in Ref. [9] is analogous to the 
Hartree approximation, and the leading order large N and Hartree approximations differ by this 
factor of 3 [1]. <j>(t) is absent in the equation of motion for Uk as the mean field is absent in Ref. 
[9], and temperature effects are ignored. 

3 Curved spacetime terms are absent in Ref. [5]. 
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dependent two-point function by solving for the mode functions Uk and then substitute them 
in the expression for the two-point function Eq. (11). Note also when comparing the above 
approaches that when one uses the large N approximation one is expanding in powers of X/N 
while in the CJT method one is expanding in powers of h. The reader is also referred to Ref. 
[14] in which the time-dependent two-point function is obtained for a system undergoing a 
phase transition with a quench albeit in Minkowski space time and without a time-dependent 
effective mass after the quench. 
Initial conditions: 

As discussed in Ref. [2], for 0(0) <C H and 0(0) ~ the growth of does not contribute 
much to the evolution of the inflaton. We assume 0(0) =0 and 0(0) = which implies that 
the mean field does not evolve with time. This is consistent with the argument presented in 
Ref. [11] that for a symmetric potential in de Sitter space the expectation values of the fields 
will be in the absence of sources because the finiteness of the effective spatial volume due 
to the presence of an event horizon allows for tunneling between the states with expectation 
value and —0. 4 Note that if (pit) = then the field <3> can be represented as an iV-plet of 

N 

„' " > 

the form ip (1, 1, ...). However, in the large N limit, working with an iV — 1-plet as we do 
gives the same results. 

We assume that the initial size of the inflating region is 0(H~ V ) and only include modes 
with k/a > 2H(0) in our expressions for (-0 2 ), energy e, pressure p, etc. [16]. (In Refs. [1,2] 
all modes are included. Our numerical results indicate that there is no significant effect of 
our lower momentum cutoff.) Then, as discussed further below, the initial conditions for Uk 
for modes with k/a > 2H(0) are given by 

4 Timescales for this tunneling and thermal initial states are not discussed in Ref. [11]. The non- 
zero overlap between the vacua associated with the two minima of our potential when the volume 
is finite may also be deduced from the discussion in Ref. [15]. It is because the vacua are related 
by a non-zero transformation when the volume is finite. 
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U k (0) = -== ; U k (0) = [-iW k0 - H(0)} U k (0), (14) 

where 

Wl, = k 2 + aSK-l + r 2 )m 2 + (£ - ±)tt(0) + ^0 2 (O) + ^ 2 (0)>] . (15) 

The term r 2 m 2 represents the corrections to the initial inflaton mass, predominantly thermal 
corrections, due to the interactions of the inflaton with other fields. We assume |(^ 2 (0)) is 
small compared to other terms in W^ and so ignore it while entering the initial conditions in 
our numerical programme. Our choice of the factor 2 in the lower momentum cutoff ensures 
that the initial frequencies defined above are real with our choice later of H(0) = 2m for all 
relevant values of k irrespective of the value of r, and that the contributions related to the 
initial fluctuations in ip can indeed be ignored in the initial frequencies. (The existence of 
the low momentum cutoff allows us to avoid the issue of imaginary frequencies at t — for 
lower momentum modes, which has to be dealt with in Ref. [2].) 

Refs. [1,2] differ in the initial conditions for the mode functions, as discussed in the 
section on the WKB approximation. Refs. [1] and [2] also differ in the initial temperature 
taken in the initial conditions for the inflaton. In Ref. [1], the authors studied a quenched 
approximation from an initial temperature T > T c to T ~ 0. The initial conditions for the 
modes of the inflaton field are specified at temperature T ~ 10 7 m R . In Ref. [2], the initial 
temperature is taken to be 0. As was noted in Sec. IVA of Ref. [1], this gives a much longer 
duration of inflation. T enters in both the coth function and the choice of r. Below we show 
that it is the influence via the former that causes a change in the duration of inflation. This 
may also be deduced from the arguments of Ref. [1]. 

III. THE LARGE K SOLUTION AND RENORMALISATION 

WKB: 

We now discuss the WKB approximation used to obtain the mode functions and their 
derivatives for large values of k. The renormalised (ip 2 ), e and e + p are then obtained by 
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subtracting the contribution of the large momentum modes. Our calculation is similar to 
that of Ref. [17]. We provide a complete derivation so as to include some details not provided 
in the literature. Furthermore our results differ from those in the literature and we point 
out the differences. 

The large k mode functions for a scalar field $ with a quartic potential were obtained in 
Ref. [7] for evolution in Minkowski space and later in Refs. [6] and [1] for evolution in FRW 
and de Sitter spacetimes respectively. These were reobtained in Refs. [17,2] with a different 
set of initial conditions than in Refs. [6,1]. The new initial conditions for the mode functions 
were obtained by first defining them for the field modes in conformal spacetime. With the 
new initial conditions, Ref. [17] focussed on radiation and matter dominated backgrounds 
while the inflationary era was studied in Ref. [2]. Below we obtain the large k mode functions 
for $ with initial conditions similar to those in Refs. [2] and further generalise the initial 
conditions to include high temperature corrections to the initial effective mass, i.e., the 
r 2 m 2 term in W k0 . The latter was included in Ref. [1] but not in Ref. [2]. In Refs. [17,2] 
it was stated that the large k behaviour of the mode functions, which is used to define the 
counterterms for renormalisation, is independent of initial conditions if one chooses the new 
initial conditions. As we see below, this is true only if temperature corrections to the inflaton 
mass at the onset of inflation are 0. 

To obtain the large k mode functions, we work in conformal spacetime. Defining the 
conformal time r from dr = dt/a(t) the mode functions in conformal time are defined via 
Uk{t) = fk( T )/C( T ) where C{t) = a(t) when r and t are related as above. We have chosen 
the normalization of the scale factor such that a(0) = C(r ) = 1. 

The equation of motion for f k is (see the Appendix of Ref. [17] for more details) 

fk(r) = (16) 
with 

M 2 = C*[-m 2 + - \)Tl\ + ^ + ^) (17) 



— 2 +k* + M(T) 
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where the tilde refers to the conformally rescaled field, i.e., 4>(x,t) = (J)(x,t)/C(t), etc. 
R(t) = QC"(t)/C 3 (t) where ' refers to differentiation with respect to r. The above expres- 
sion for the effective mass differs from that in Eq. (A7) in Ref. [17]. 

The initial conditions for the conformal time mode functions f k are chosen to correspond 
to an initial density matrix that represents a system in local thermodynamic equilibrium and 
which commutes with the Hamiltonian in conformal spacetime at r = r [17]. Therefore, for 
modes with k/a > 2H(0), 

f k (r ) = -±= ■ f k (r ) = -iW m f k (T ) ■ W k0 = ^/k 2 + M 2 . (18) 

M.q is the initial effective mass in terms of conformal time. 

Ml = C 2 (r )[(-1 + r 2 )m 2 + (£ - \)K{t q )\ + ^ 2 (r ) + ^ 2 (r )> (19) 

Note that the temperature corrections r 2 m 2 above is not included in Eq. (17) which is 
valid at times after the quench. The presumption here is that the other fields to which the 
inflaton couples go out of equilibrium during the quench at a temperature much below Tj 
at which stage their contribution to the effective mass of the inflaton is small and can be 
ignored, and that any subsequent growth in their fluctuations during inflation does not lead 
to a significant contribution to the effective inflaton mass. 

These initial conditions for modes with k/a > 2H(0) translate to 

U k (0) = -±= ■ U k (0) = [-iW k0 - H(0)\ U k (0), (20) 

V W k0 

where W k o is W k o expressed in comoving coordinates. 

Dividing the equation of motion for f k Eq. (16) by k 2 /m 2 , we treat m/k = 5 as a small 
parameter and apply WKB theory to solve for f k . We shall seek solutions for k > 10m. 
Using a WKB solution for f k of the form [18] 

D k (r) ~exp[-=S(r)] = exp[- ^ 5^„] = e K R ^ )dr , (21) 

we obtain two solutions, D k i(r) and its complex conjugate, which agree with the expressions 
in Eqs. (A14) and (A15) of Ref. [17], namely, 
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t-> / \ -l , o f \ - R lA T ) , R 2,k( T ) - R 3,k( T ) , R ±A T ) , ^oo^ 

i4(r) = -ifc + i2o, fc (r) - — + — ~ 2 t— p— + — ^— + • • • (22) 

and its complex conjugate. The coefficients are given by 

R o,k = ; R lik = \m\t) ■ R 2 , k = - l -R' lk = ~M 2, (t) 



R s, k = \ - R l, k ) = -\m*'{t) - \m\t) 

R *,k = ~\ (K, k + 2R 1>k R 2>k ) = ±-M 2 "\t) + \m"{t) . (23) 



2 V L ' K J 8 8 

2 (-^,fc + 2R l,k R 2, k ) = — 

Imposing the initial condition for f k (ro) implies 



fk(r) = ^= [(1 + -y)D kl (r) + (1 - 7)^i(r)] • (24) 

7 is then obtained by imposing the initial condition on f k (r ). We find that the real and 
imaginary parts of 7 are given by 

7 * = 1 + P + F 

where 

91 = \{Ml-Ml + ) 

92 = \m™ -\{M%- Ml + ) - \Ml + {Ml - Ml) 

93 = Mb 

9, = ~M™ - \Ml + Ml' + + \MlMl' + (26) 

and Ml + = M 2 {tq). Note that Ml — Mq + = r 2 m 2 since Ml and Ml + are the effective 
masses in conformal spacetime just before and after the onset of inflation. 
Expressing D k i{r) as e x+lY and using the above relations we get 

\f k \ 2 = _L e 2X [cos 2 y + 7^sin 2 F-27 / cosFsinF + 7f sin 2 F] (27) 

Ultimately the above expression will be integrated over all k. Expressing trignomet- 
ric functions as exponentials there are terms in the integrand that for large k go as 

11 



k 2 k 1 l exp[±2ik(r — r ), I > 0. For large k(r — r ) we set their integrals to 0. So we 
replace cos 2 F and sin 2 Y by | and cos F sin F by and write the effective \fk\ 2 as 
1 



2W k0 
1 



k 



1 



1 M 2 {t) 1 f A^ 2 "(r) 3.M 4 (r) (Aj 

r~ + F I 8 + : 
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2 _ KAl \1 

yvi o+J 



(28) 



where we have also ignored the term proportional to since it is of order 1/k 6 and can be 
ignored in the large k limit. \f' k \ 2 is given by 
1 



\fl\ 2 = 



,2X 



W fc0 



{i? fcjR cosF - R kI sinF - ^i(R kR sinF + R kI cosF)}' 
+ ln{RkR sin F + i? fc/ cos F} 2 



(29) 



where R k R and are the real and imaginary parts of R k defined in Eq. (21). Again, the 
effective \f' k \ 2 relevant for subsequent integration over all k and excluding terms of higher 
order than 1/k 3 is 
1 



l\2 



2W, 



fcO 



= k 



1 + 



e 2X [l+l R ]R 2 kI 
1 A^ 2 (r) 



_ MM + (M 2 -M 2 + ) 2 + MHr) _ 



(30) 



k 4 I 8 8 8 4 

(We obtain \ f k \ 2 to 0(1/ k 5 ) and \f k \ 2 to 0(l//c 3 ) since the former is multiplied by a factor 
of k 2 in the expression for energy and pressure. However in the expression for (ip 2 ) we only 
keep terms to 0(l/k 3 ) in \f k \ 2 .) 

We reiterate that the terms in Eqs. (27) and (29) set to 1/2 and attain these values 
only after integration over all k. There is no mention of this in the literature. Furthermore, 
during inflation r = —l/(aH). Treating H as a constant equal to 2m we can easily check 
that for k > 3m, k(r — r ) >> 1 is valid within two e-foldings or so of inflation. Once 
inflation ends the universe expands as t n (n < 1) and r = j^i^aH)^ 1 . Then r is positive 
and increases with time. Hence k{r — t ) >> 1 will continue to be valid. However for 
k < 2m, k(r — ro) >> 1 will never hold during the inflationary era. Therefore later we shall 
use the above large k expressions only for k > 10m R , where m R is the renormalised mass 
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parameter, rather than for k > rriR as in Refs. [1,2], while subtracting off the divergences 
during renormalisation. 

As mentioned earlier, one can see that it is only for r = 0, i.e. when = J\A% + , that 
the large k subtractions needed to renormalise the divergences are independent of initial 
conditions. For r = 0, our result for |/fc(T)| 2 agrees with Eq. (A24) of Ref. [17]. Our result 
for \ f' k {r)\ 2 disagrees with Eq. (A25) of Ref. [17]. However it agrees with the corresponding 
expression in Eq. (5.24) of the review article [19] indicating a possible typographical error 
in Ref. [17]. (Note that the expression for \fk{r)\ 2 (r = 0) in Eq. (5.24) of Ref. [19] disagrees 
with both our expression and that of Ref. [17] indicating a possible typographical error in 
Ref. [19].) 

The corresponding expressions for the comoving time mode functions are 
\h(r)\ 2 



\U k {t)[ 



C 2 (r) 
1 



(31) 



ka 2 (t) 2k 3 a 2 (t) 



B(t) 



\Uk(t)\' 



8k 5 a 2 (t) 
S (2) + 0(l/k 7 ) 

M' k {r)\ 2 
C 2 (r) 
1 



3B(t) 2 + a(t)a(t)B(t) + a{t) 2 B{t) + (B - B 0+ f 



+ 0(l/k 7 ) 



1 

CHr) 
k 

a A (t) + 2ka 4 (t) 
1 

4 



+ [ H — 



H d s 



\fk{r) 



(32) 
(33) 



B(t) + 2a(tf 



8k 3 a 4 (t) 
+2B(t)(B -B 0+ ) 

S {1) + 0(l/k 5 ) , 



B(t) 2 - a(t) 2 B(t) + 3a(t)a(t)B(t) - 4a 2 (t)B(t) + (B - B 0+ ) 



+ 0(l/k 5 ) 



(34) 



where B(t) is M. 2 {r) expressed in comoving coordinates and B and B 0+ are similarly 
defined. 

Renormalisation: 

As in Ref. [1,2] we use the above expressions to write the renormalised (ip 2 ), (^ 2 )_r, as 
(hereafter the subscript R refers to renormalised quantities) 
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(35) 

k represents a renormalisation scale. As mentioned above, the large k expressions used in 
the subtraction are valid for k > lOm R (after two e-foldings). Hence we take k to be 10m R 
rather than m R as in Ref. [2]. We also choose to not rescale (ip 2 ) R to be at t = 0. As in 
Ref. [1] we shall take H(0) = 2m R . 

B(t) is now expressed in terms of renormalised parameters making use of the condition, 
M R (t) = M 2 (t) discussed in Sec. IV of Ref. [6]. (Note that M 2 (t) is written in terms of 
bare parameters.) Therefore 

B(t) = a 2 (t)[-m% + - \)K{t) + ^0 2 (t) + ^ 2 (t)) R ] (36) 

Since B contains (^> 2 )r one has to solve Eq. (35) to obtain an explicit expression for (tp 2 ) R . 

Besides the presence of the first 6 function, our approach above also differs from that 
of Ref. [1,2] in that the second function applies to all terms in the subtraction and not 
just to the last term that suffers from an infrared singularity. 5 The difference between the 
subtractions under the prescription above and that of Refs. [1,2] is finite and so both are 
valid. However, as we discuss below, a prescription that subtracts all terms in an expansion 
for \Uk\ 2 for large k modes makes it simpler to obtain the expressions for the renormalised 
energy and energy plus pressure that satisfy the energy conservation equation. 

Similar to the approach in Refs. [1,2], the renormalised energy and energy plus pressure 



5 The expression for (ip 2 (t)) R in Eq. (3.18) of Ref. [1] is missing one subtraction term (and the 
coupling). We thank D. Cormier for confirming that this is a typographical error. There are similar 
errors in Eq. (3.20) of Ref. [1] and in the expression for the renormalised fluctuations in Sec. V of 
Ref. [2]. 



/rl°k 
-^coth[H/ fc0 /(2T l )] 
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are defined as 



6 7 




1 r H s k r 
+ 2 J 2^coth[W fe o/(2T,)] [\U k \ 2 Q(k - 2H(0)) - S^Q(k - k) 

W kR {\U k \ 2 Q{k - 2H(0)) - 5< 2 >e(* - «))] - 



(37) 



N 



2 + / 2^cotfr[WV(ZTi)] [|t/ fc | 2 9(fc - 2ff(0)) - 5«0(A; - «) 
+ ^(l^| 2 ©(^ - 2//(0)) - 5( 2 )0(A; - «)) 



(38) 



The renormalised energy and pressure should satisfy the energy conservation equation 



For the energy conservation equation to hold the 6 function for the subtracted terms 
should be chosen compatible with the subtraction in {^ 2 )r. This has been done in Ref. 
[1] and above. Refs. [1,2] differ from each other in their subtractions for the energy and 
energy plus pressure and we find that the expressions for and (p + e)r in Ref. [2] do 
not satisfy the energy conservation equation because the lower momentum cutoff for the 
different subtracted terms are not chosen appropriately. We find that a simple prescription 



6 There is an error in the last term in the expression for the renormalised energy density in Eq. 



7 An alternate prescription of obtaining the renormalised energy-momentum tensor using adiabatic 
subtraction is provided in Ref. [20], albeit for chaotic inflation. In this prescription one expands 
bare quantities in powers of the derivatives of lna which is equivalent to expanding about solutions 
for Minkowski spacetime. 



e R (t) + 3h(t) (p + e)«(t) = 0. 



(39) 



(3.21) of Ref. [1]. 
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is to adopt the same lower momentum cutoff for all the subtracted terms in {ip 2 )n, €r and 
(p + e) R , and not to rescale (^> 2 )r at t = to be 0. 

An issue that appears to have been overlooked in Refs. [17,2] is that the WKB solutions 
obtained above are also valid only for k 2 3> M 2 . In fact, one can see that if one guesses a 
solution for Eq. (16) of the form exp[i J(k 2 + M 2 )^ dr] then the exponent will approximate 
to a series of the form Eq. (22) (in the limit of slowly varying A4 2 ) only if M 2 <C k 2 . 
Because of the factor C(r) in M this condition breaks down very quickly for the lower 
momentum modes in the integrals over k. Alternatively, the validity condition for the WKB 
solutions is written as 5 n S n+ i <C 5 n ~ 1 S n as 5 — > [18], i.e. succeeding terms in the series 
in the exponent of the WKB solution should be smaller than previous ones. Applying this 
condition to the non-zero terms in the series for R k we find that it holds only for k 2 3> M 2 . 
Similarly, in the expansion for \Uk\ 2 in powers of \jk in Eq. (3.2) of Ref. [1], the second 
term is smaller than the first only if k 2 /a 2 3> M 2 (t) — 1Z/6, and similarly for the second 
and third terms. Thus the subtractions being carried out do not match the true solution 
for increasingly larger numbers of modes as inflation progresses. (Though the WKB series 
will not approximate the true solution for k < M. this may not interfere with the energy 
conservation equation as the individual terms in the series are obtained consistently order 
by order from the equation of motion.) 

Now, for the inflaton field the dominant contribution to (ip 2 )n, €r and pr come from low 
momentum modes. When a mode becomes unstable, i.e., when its effective frequency u\ R 
becomes negative, it grows exponentially fast. Obviously low momentum modes become 
unstable earlier and hence start growing earlier. Furthermore, after a few e-foldings of 
inflation the contribution of additional modes that become unstable does not add much 
to the contribution of the modes that left in the first few e-foldings. As discussed above, 
for the dominant low momentum modes the WKB solutions are not valid approximations 
of their behaviour. However, as seen in Fig. 2, the subtractions carried out in Refs. [1,2] 
based on the WKB solutions do not affect the contribution of these modes much. This may 
explain why the numerical solutions of Ref. [2] did not show any conflict with the energy 
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conservation equation even though the counterterms were subtracted improperly. Note that 
for a field with a positive mass squared, as also considered in Ref. [17], this explanation will 
not hold. 

A comparison with the standard prescription 

We now compare the renormalisation schemes adopted in Refs. [1,2] based on Ref. [6] 
with the approach outlined in Ref. [9]. (See also Ref. [16].) The basic difference is that in 
Ref. [9] all modes that have left the horizon during inflation are included in the integral over 
k in (V ;2 )- This provides an ultraviolet cutoff of Ha for the integral. (In this sub-section 
we assume H is constant during inflation for easy comparison with Refs. [9,16] where this 
is also assumed.) In the scheme adopted in Refs. [1,2] all divergent terms in the WKB 
solution obtained above are subtracted for modes k> k from \U k \ 2 to obtain (ip 2 )n. In the 
former scheme more and more modes contribute to the integral in (iP 2 )r as they cross the 
horizon but this does not seem to be the goal in the latter scheme. But it is precisely the 
contribution of the newer modes leaving that gives the standard time dependence of {ip 2 )n 
for a massless or massive field, as in Ref. [21]. 

For example, in the case of a massless field for which we know the exact solution for £4, 
the prescription of Ref. [9] gives 8 



(The second term is ignored as it is suppressed by a 2 .) Presuming that the intent in the 
scheme of Ref. [6] was to subtract the contribution of modes greater than some renormali- 
sation scale k, and making use of the solution obtained in Ref. [9] to correctly subtract the 
contribution of all modes k > k, {4> 2 )r will be 



8 As stated in Ref. [9], this approach provides a simple but intuitive way of obtaining (iP 2 )r. A 
more rigourous derivation with proper regularisation and renormalisation is given in Ref. [22]. 




HH 



(40) 



4tt 2 ' 
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(2it) 3 Jh 2k 3K H 2 a 2 ' 

« §^Hk/H) . (41) 

(We have chosen the lower momentum cutoff to be H in keeping with Ref. [9].) One sees that 
the two prescriptions give very different results which must be kept in mind while calculating 
the effect of fluctuations. 

For our field ip with a negative mass squared the renormalisation prescription based on 
Ref. [6] implies 

. , 9 . f°° 47rk 2 dk\U k (t)\ 2 , rTTr _ X1 f°° 4irk 2 dk\U K (t)\ 2 , rTTT ., _ X1 . 
= ( 27 r)3 2 coth[^o/(2T,)]-/ K (27r) a 1 coth[^ /(2TO] (42) 

« 4vrA; 2 dA;|[4(t)| 2 



2// (2tt) 3 



-coth[W fc0 /(2T i )] 



+ / k 2 (|f/fcW|2 " I^WI 2 )coth[^ /(2T l )], (43) 

where Uxit) are the WKB solutions and |t/x(<0| 2 above contains terms of order 1/k and 
1/k 3 . As discussed earlier, the WKB solutions U K (t) in the second integrals are not an 
approximation for the corresponding Uk(t) for k < \M\. (Eq. (43) is written for times when 
K < which is all times if k = ttir and after 2 e-foldings if k = lOm^.) On the other 
hand, the renormalisation scheme of Ref. [9] would imply 

. f Ha 4irk 2 dk \U k (t)\ 2 , rTTr 

^ * = / fo ^ o coth Wiw/ 2T, (44) 

(for times when _£/a > 2H, i.e., after 1 e- folding). 

Both expressions for (ip 2 )n appear different. (Since they differ by a finite amount, in 
principle, both are valid.) However, in Fig. 1 we have plotted {ip 2 )n obtained under both 
renormalisation schemes and they are similar. For the case of the inflaton with a negative 
mass squared lower momentum modes are growing modes and the contribution of the lowest 
few modes k < 20m# dominates the integrals in (iP 2 )r. For these modes the WKB sub- 
tractions do not affect their contributions, as seen in Fig. 2. Therefore once \A4\ and Ha 
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become larger than 20m R , which happens within 3 e-foldings, (ip 2 )n obtained under both 
renormalisation schemes are similar. 
Some additional comments: 

In Ref. [1] the authors cite the possibility using the growth of fluctuations of the inflaton 
as a new means of evolving the inflaton in its potential and providing a graceful exit from 
the inflationary phase. However, following the evolution of the fluctuations rather than 
the mean field to study the inflationary phase transition has been studied earlier in Refs. 



In Sec. V(D) of Ref. [2] the authors introduce the notion of a zero mode assembly 
consisting of super-horizon sized low momentum modes, which acts as an effective zero 
mode after a few e-foldings of inflation even in the absence of <f>(t). The authors have also 
discussed the classicalisation of these low momentum fluctuations. We point out that this 
is similar to ideas discussed in Sec. IIIA of Ref. [23] and in Sec. 8.3 of Ref. [16] in which 
the early evolution of the inflaton is due to the growth of quantum fluctuations (i/) 2 )r but 
after some time t c i one can study inflation by solving the classical equation of a motion for a 



purely time dependent field with the initial value at t d taken (in Ref. [16]) to be (^> 2 (t d )) R . 

IV. NUMERICAL RESULTS 

In this section we present our results after numerically solving the large set of differential 
equations representing a continuum of ^-values. The numerical code is based on fourth order 
Runge-Kutta for evolution of the mode functions and IMSL algorithms for interpolation of 
discrete data and integration to obtain (iP 2 )r. We have chosen A = 10~ 12 consistent with the 
calculation of density perturbations in Ref. [2] and we set = 0. For numerical purposes 
we follow Ref. [1] in introducing the following dimensionless quantities 



[21,9,11]. 



t = m R t 



h 



m R 



H 



k 

Q = 



(45) 



Tu = Ti/ m R . 



(46) 
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8tt 2 
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We treat H as a constant. Assuming m R ~ 10 13 GeV we set h = 2. 

In Fig. 1 we present the growth of gE under the renormalisation schemes of Ref. [6] and 
Ref . [9] . The plot following the renormalisation scheme of Ref. [6] contains the modifications 
suggested above, namely, only modes k/a > 2H(0) are included in the integrals over k, the 
theta function for the subtracted quantities applies to all terms being subtracted and the 
renormalisation scale k is taken to be 10m#. 9 It is not possible to distinguish between the 
two curves and the quantum fluctuations grow dramatically between i = 60 and 100 and 
approach an asymptotic value of 1. 

In Fig. 2 we plot gT, with and without the subtraction of Eq. (42). The two curves lie 
on top of each other. This implies that the effect of the subtraction is negligible. 

In Fig. 3 we plot g£ for different values of T u . The initial temperature Tj enters in 
two places in the calculation of the inflaton dynamics in Ref. [1]. Firstly, it appears in the 
temperature dependent correction to the effective mass for ip in the initial conditions (r is 
taken to be Tj/T c in Ref. [1]) and secondly, it appears in the expressions for (V> 2 ), energy, etc. 
in the factor coth[Wfco/ (27*)] ■ To distinguish between these two effects we choose Tu = 0, 200 
and 10 7 and take r = for the first case and r = 2 for the other cases. For T u = 0, g£ 
begins to grow appreciably at i « 70 and achieves an asymptotic value of 1. This agrees well 
with the results in Ref. [2]. For Tu = 10 7 , gT, begins to grow appreciably at i « 30. This 
again agrees well with the results in Ref. [1]. The corresponding timescale for T u = 200 is 
60. Comparing the plots for Tu = 200 and 10 7 with the same value of r we see that it is the 
presence of the coth term representing initial thermal fluctuations of the i/j field that plays a 
substantial role in modifying the duration of inflation. As a further test we have also varied 
r in Who till r = 10 and seen that there is no significant dependence on r. 

The inflationary phase transition is completed when the effective mass for the inflaton 
modes — 1 + gT, achieves a value of 0. On the other hand, the inflationary era ends when 



Subsequent figures follow the renormalisation scheme of Ref. [6] with the same modifications. 
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a < 0. Now a = — -^(e^ + 3p#) and hence inflation ends when p R becomes greater than 
— \zr- Since H / H 2 = — 1(1 + 2s), inflation ends when H / H 2 becomes less than -1. Since 
we do not treat the Hubble constant dynamically we take the breakdown of the constant 
H assumption, when the fluctuations begin to grow dramatically, to be an indicator of the 
end of inflation. As the fluctuations grow the universe will transit from an exponentially 
inflating universe to a power law inflating universe and will then enter the reheating era. 
Therefore, hi e , where t e is the time when the fluctuations begin to grow dramatically, gives 
us a lower bound on the number of e-foldings of inflation. 

Understandably, the existence of frozen-in initial thermal fluctuations in ip decreases the 
time required for the total fluctuations (ip 2 )n to reach its final value. If one assumes that the 
inflaton was in thermal equilibrium prior to the onset of inflation then it is likely it will go out 
of equilibrium soon after the onset of inflation and the above shows that it is important to 
include the thermal fluctuations of the initial state while estimating the duration of inflation. 
If the inflaton goes out of thermal equilibrium earlier than the onset of inflation then the 
contribution of the thermal fluctuations frozen in at the higher temperature can be even 
larger. 

V. CONCLUSION 

In this paper we have firstly compared the approach of Refs. [1,2] with that of earlier 
works to study the dynamics of the inflaton. We point out that for the Hartree and the lowest 
order large N approximations employed in Refs. [1,2] with the CTP formalism the equations 
of motion are akin to those obtained in standard field theory. We have reobtained the large 
momentum mode solutions for the inflaton field using the WKB method and pointed out 
some differences with earlier results. Because of these differences and the very involved 
nature of these calculations we have included many explicit details of the calculations. We 
have also included temperature corrections in the initial effective mass of the inflaton field 
and pointed out that the assertion in Refs. [17,2] that their choice of initial conditions 
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in conformal spacetime gives counterterms that are initial conditions independent is only 
true if there are no thermal corrections to the initial inflaton mass. We have discussed 
different aspects of the renormalisation prescription adopted in Refs. [1,2]. In particular, we 
have pointed out that the asymptotic expressions used in Refs. [1,2] for the large momentum 
modes during renormalisation are valid for large k/a and not large k. We have also compared 
the renormalisation prescription of Refs. [1,2] with that of earlier works and pointed out 
their differences. Despite these differences the results of these different approaches are the 
same because for an inflaton with a negative mass squared it is the low momentum modes, 
and not the high momentum modes regularised by renormalisation, that give the dominant 
contribution. 
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Figure Captions: 

Fig. 1: Renormalised fluctuations of the inflaton field, g£, are shown as a function of 
normalised time t for the two renormalisation schemes discussed in the text (with constant 
H). r = 2, h = H/m R = 2 and the normalised initial temperature T„ = Ti/m R = 200. The 
solid line corresponds to the (modified) renormalisation scheme of Ref. [6] while the points 
correspond to the renormalisation scheme of Ref. [9] . 

Fig. 2: Fluctuations of the inflaton field, with and without the WKB subtractions (solid 
line and circles respectively), are shown as a function of normalised time t. All parameters 
are as in Fig. I. 

Fig. 3: gY, vs. t for T u = 0, 200 and 10 7 . r = 2 for T u = 200 and 10 7 and r = for 
Ta = 0. h = 2 for all plots. 
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Figure 1 
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